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Abstract 

The fractional diffusion equation is derived from the master equation of continuous 
time random walks (CTRWs) via a straightforward application of the Gnedenko- 
Kolmogorov limit theorem. The Cauchy problem for the fractional diffusion equa- 
tion is solved in various important and general cases. The meaning of the proper 
diffusion hmit for CTRWs is discussed. 



1. INTRODUCTION 

This paper provides a short, but self-contained, introduction to fractional diffusion. The 
readers will find the basic ideas behind the derivation of the fractional diffusion equation 
starting from continuous-time random walks. We have included formulae for the solution 
of the Gauchy problem which can be numerically implemented and used for applications. 
Special care has been used to avoid unessential mathematical technicalities. Even if far 
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from exhaustive, the bibhography should give a sufficient number of entry points for further 
reading. The following sections are based on a series of papers about the application of 
fractional calculus to financeB'iHS. 

It is our hope that theoretical and experimental condensed matter physicists will find 
this work useful. 

The paper is divided as follows. In Section 2, we outline the theory leading to the 
time- fractional master equation. In Section 3, the transition to the space-time fractional 
diffusion equation is discussed. Section 4 is devoted to the solutions of the Cauchy problem 
for the fractional diff'usion equation. The main results are briefly summarized and discussed 
in Section 5. In Appendix A and Appendix B, we introduce the definitions of fractional 
derivatives in time and space, respectively, entering the fractional diffusion equation. 



2. STEP ONE: TRANSITION TO THE TIME-FRACTIONAL MAS- 
TER EQUATION 

Let X be the position of a diffusing particle in one dimension. Let us assume that both 
jumps = x{t.i) — x{ti-i) and waiting times between two consecutive jumps ti = ti — U-i 
are i.i.d. random variables described by two probability density functions: w{^) and iP{t). 
According to the model of continuous-time random- walk (CTRW), introduced by Montroll 
and WeisBi, the evolution equation for p{x, t), the probability of finding the random walker 
at position x at time instant t , can be written as follows, assuming the initial condition 
p{x, 0) = 6{x) (i.e. the walker is initially at the origin x = 0)i, 



p{x,t) = (5(x)^'(t)+ / Tpit-t') 
Jo 



w{x — x') p{x' , t') dx' 



dt' . 



(2.1) 



where 

v(t') dt' = i- i^it') dt' , m = --^{t) . (2.2) 

The master equation of the CTRW can be also derived in the Fourier-Laplace domain. 

The integral Jq tp{t') dt' represents the probability that at least one step is taken at some 
instant in the interval [0,r). Thus, ^'(t) is the probability that the diffusing quantity x 
does not change value during the time interval of duration r after a jump. 

In a paper by Mainardi et alJl, an alternative form of Eq. (2.2) was presented in terms 
of a convolution between the first time derivative of p{x, t) and a suitable kernel. The 
resulting equation can be interpreted as an evolution equation of generalized Fokker-Planck- 
Kolmogorov kind. It reads: 

rt Q r+oo 

/ <^(t -t') —p(x,t')dt' = -p(x,t) + w(x - x')p(x',t)dx' , (2.3) 

Jo of J-oo 

where the "auxihary" function $(t) is such that ^{t) = /g $(t - t') ip{t') dt' . Eq. (2.3) can 
be obtained by Fourier-Laplace transforming Eq. (2.1) and by suitable assumptions on the 
Laplace transform of the function $(t). 

In general, a CTRW is a non-Markovian process. A CTRW becomes Markovian if (and 
only if) the above memory function is proportional to a delta function so that ^{t) and 
V'(i) differ only by a multiplying positive constant. By an appropriate choice of the unit of 
time, we can write $(i) = 6(t) , t > . In this case, Eq. (2.3) becomes: 
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/ + 00 
w{x — x') p{x' ,t) dx' , p{x,0) = 6{x). (2.4) 



d_ 

Up to a change of the unit of time, this is the most general master equation for a Markovian 
CTRW; Saichev & Zaslavsky cah it the Kolmogorov-Feller equatio'M. 

Eq. (2.3) allows a natural characterization of a peculiar class of non-Markovian processes, 
where the memory function, <I>(t) has power-law time decay. Within this class, an interesting 
choice is the following: 

In this case, $(t) is a weakly singular function that, in the limit /?—>!, reduces to 
$(t) = 6{t) , according_to the formal representation of the Dirac generalized function, 
6{t) = t~^/T{0) , t > (ji. As a consequence of the choice (2.5)1, Eq. (2.3) can be writ- 
ten as: 



-—-np(x,t) = —p(x,t)+ / w(x — x') p(x' ,t) dx' , p{x,0) = 6(x) , (2.6) 

OtP J_oo 

where /dt^ is the pseudo-differential operator explicitly defined in the Appendix A, that 
we usually call the Caputo fractional derivative of order [3 . Eq. (2.6) is a time- fractional 
generalization of Eq. (2.4) and can be called time- fractional Kolmogorov-Feller equation. 

Our choice for <I>(t) implies peculiar forms for the functions ^{t) and %l){t) generalizing 
the exponential behaviour of the waiting time density in the Markovian case. In fact, we 
have for t > oi: 

^{t)=Ep{-tP), ^(t) = -jEp{-t^), 0</3<l, (2.7) 

where Ejj denotes an entire transcendental function, known as the Mittag-Leffler function 
of order (3 , defined in the complex plane by the power series 

oo 72 

^.M-Ef,^. '3>0. ^^C. (2.8) 

Detailed information on the Mittag-Leffler-type functions is available in the literatureilli0. 

Erom the properties of the Mittag-Leffler function, it can be shown that the corresponding 
survival probability and waiting-time pdf interpolate between a stretched exponential, for 
small waiting times, and a power-law decay, for large waiting times. Such a behaviour has 
been observed in Mainardi et al.i and in Raberto et al.i, where the function ^'(t) has been 
estimated from empirical financial data. 

As a final remark, it is important to notice that different choices of the kernel ^{t) in 
eq. (2.3) are possible, leading to different properties of the waiting-time pdf and different 
generalized Kolmogorov-Eeller evolution equations for p(x,t). 
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3. STEP TWO: THE DIFFUSION LIMIT 

In the physical hterature, many authors have discussed the connection between continuous- 
time random walks and diffusion equations of fractional order, with different degrees of 
detaiii-§. A sound proof of the equivalence between fractional diffusion and CTRW has 
been given by Hilfer & Antonlll. However, in order to perform the transition to the diffusion 
limit, we shall use a different approach. We shall start from Eqs. (2.4) and (2.6), and 
pass through their Fourier-Laplace counterparts. The stochastic process whose probability 
density evolves according to those equations is a random walk originating from a sequence 
of jumps, each jump being a sample of a real random variable Y. During the time interval 
tn ^ t < tn+i , the particle position is Yi + y2 + ■ ■ ■ ■ The Yk are i.i.d. random variables 
all described, as Y, by the pdf w{x) . Let us denote by w{k) the characteristic function 
corresponding to the probability density w^x). 

Let us specify some conditions on the pdf w{x). The requirement is that, if a = 2 : 

/ + 00 

x'^ w{x) dx < oo , (3.1) 

whereas, if < a < 2 : 

w{x) = {b + e{\x\)) \x\-^"'^'^\ b>0, e(|x|) ^ as ^ oo . (3.2) 

In Eq. (3.2), 6 > and e{\x\) is bounded and 0(|j;|~'') with r/ > as — > oo . Let us 
furthermore recall the necessary requirements w{x) > 0, and the normalization condition 
J^^ wix) dx = l. 

Let us now consider a sequence of random process pdf's ph{x,t) describing scaled jumps 
of size h Yk instead of Yk , with a speed increase of the process by a factor {the scaling factor) 
ii^^'/P ^ where ^ must satisfy some conditions which will be specified later. The pdf 
of the jump size is Wh{x) = w{x/h)/h , so that its characteristic function is Wh{n) = w^nh) . 
For < a < 2 and < /3 < 1 , Eq. (2.6) (including (2.4) in the special case /? = 1) is 
replaced by the sequence of equations 

QP r+oo 

l^h" -^Ph{x,t) = -ph{x,t) + J Wh{x - x')ph{x',t)dx' . (3.3) 

By Fourier-Laplace transforming and by recalling the Laplace transform of the Caputo 
time-fractional derivative, defined by Eq. (A. 2), we have 

M/i" {s^j^hiK,s)-s^-'] = [whit^)-l]Phi^^,s). (3.4) 

We shall now present arguments based on the classical central limit theorem or on the 
Gnedenko limit theorem, (see the book by Gnedenko & Kolmogorov@) both expressed in 
terms of the characteristic functions. The Gnedenko limit theorem is a suitable generaliza- 
tion of the classical central limit theorem for space pdf's with infinite variance, decaying 
according to condition (3.2). 

The transition to the dijfusion limit is based on the following Lemma introduced by 
Gorenfloii: 
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With the scaling parameter 



■f 9 

— , II a = 2 , 

2 ' 

if < a < 2 , 



(3.5) 



, T{a + 1) sin(Q7r/2) 



we have the relation 



lim^^^^ 1 = -|k|°, 0<a<2, kGR. (3.6) 



Now, it is possible to set 

nui K.\ = — 

flh 



w{Kh) - 1 

Ph{>^) = rz — , (3-7) 



and the sequence of equations (3.4) reads 

s^p'h{K,s)-s^~'^=Ph{K)Ph{>^,s). (3.8) 

Then, passing to the hmit /i ^ 0, thanks to (3.6), we get: 

/^(k,s) - /"^ = -|k;|"^(«;,s) , 0<a<2, 0</3<l. (3.9) 

By inversion and using the Fourier transform of the Riesz space-fractional derivative, defined 
in Eq. (B.2), we finally obtain the equation: 

-Q^Poix, t) = -Q^Poix, t) , po{x, 0) = 5{x) , (3.10) 

which is a space-time fractional diffusion equation. In the limiting cases /3 = 1 and a = 2, 
Eq. (3.10) reduces to the standard diffusion equation. 

We have presented a formally correct transition to the diffusion limit starting from the 
general master equation of the CTRW, namely Eq. (2.1) or Eq. (2.3). By invoking the 
continuity theorem of probability theory, see e.g. the book by Lukac^^, we can see that the 
random variable whose density is ph{x,t) converges in distribution ("weakly" or "in law") 
to the random variable with density po{x, t) . 

Solving (3.8) for ph{K, s) , and (3.9) for po(k, s) , gives: 

^(^'^^ = ./3-p,(k) ' ^("'^^ = i^r^' (3-11) 

which yields: 

pi,{K,t) = Ep [phiK)t^) , m^^,t) = Ep (-i^rt^) . (3.12) 

By (3.6) Ph{i^) as /i — > , hence 

Ph{x,t) ^ PQ{x,t) , for t>0, /i^O. (3.13) 
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4. SOLUTIONS AND THEIR SCALING PROPERTIES 

For the determination of the fundamental solutions of Eq. (3.12) in the general case {0 < 
a < 2 , < (3 < 1} the reader can consult Gorenflo et al.0 and Mainardi et alS. We 
also refer to the above references for the particular cases {0 < a < 2 , f3 = 1} and {a = 2 , 
< P < 1}, already dealt with in the literature. 

For parameters in the interval < a < 2, and < /3 < 1, the Cauchy problem in Eq. 
(3.12) can be solved by means of the Fourier-Laplace transform method. 

The solution (Green function) turns out to be: 

P0i^'t) = 7LWe.,J-^)- (4.1) 



^f3/a ""./^ y^l3/. 

The function VFQ^^(ti) is the Fourier transform of a Mittag-Leffier function: 



1 f + OO 

Wo.A^) = - e-^^^Ep{-\qndq. (4.2) 



Indeed, E/^ is the Mittag-Leffier function of order P and argument z 
In the limiting case < a < 2 and /? = 1, the solution is: 



Poi^, t) = TuzLa ( ttt;: ) ' (4-3) 



where La{u) is the Levy standardized probability density function: 



27r 

whereas, in the case a = 2, < /? < lH, we get 



1 r+oo 

Lain) = — I e-*^"-l<?l dq, (4.4) 

oo 



PO 



(^>*) = ^*-"/'^/3/2(^) , (4.5) 



where denotes the M function of Wright type of order (3/2 . 

Remarkably, a composition rule holds true, and it can be shown that the Green function 
for the space-time fractional diffusion equation of order a and /? can be written in terms 
of the Green function for the space-fractional diffusion equation of order a and the Green 
function for the time-fractional diffusion equation of order 2/311: 

(x, t) = t-f^ r r-i/" La (x/r^/") (r/t^) dr . (4.6) 



Pol 

Finally, as written before, in the case q = 2, /? = 1, Eq. (3.12) reduces to the standard 
diffusion equation, and the Cauchy problem is solved by: 

Po(x, t) = _i_ exp(-xV(4t)) = t-^" G , (4.7) 

where G{x) denotes the Gaussian pdf 

G(x) = -^exp(-xV4). (4.8) 
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5. SUMMARY AND DISCUSSION 

Applications of fractional diffusion equations have been recently reviewed by Uchaikin &; 
ZolotarevS and by Metzler & KlafteiS. After that, other contributions appeared on this is- 
sue, among which we quote the paper of Zaslavsky in the book edited by Hilfeiill, the papers 
by Meerschaert er al@ and by Paradisi et al.H and the letter by West and NonnenmacheiS. 

In this paper, a scaling method has been discussed to get the transition to the diffu- 
sion limit in a correct way, starting from the CTRW master equation describing the time 
evolution of a stochastic process. Moreover, the solutions of the Cauchy problem for the 
fractional diffusion equation have been listed for the various relevant values of the fractional 
derivative orders a and (3. 

Various formulae which can be useful for applications have been presented. In principle, 
given a diffusing quantity, the waiting-time density, the jump density, and the probability of 
finding the random walker in position x at time t are all quantities which can be empirically 
determined. Therefore, many relationships presented above can be corroborated or falsified 
in specific contexts. 

As a further remark, it may be useful to add some comments on the meaning of the 
diffusion limit taken in section 3. 

The factor fih" can be viewed as causing the jump process to run faster and faster (the 
waiting times becoming shorter and shorter) as h becomes smaller and smaller. Replacing 
the density by the density Wh{x) = w{x/h)/h , and, accordingly, the jumps Y by hY , 
means that the jump size becomes smaller and smaller as the scaling length h tends to zero. 

An alternative interpretation is that we look at the same process with a discrete number 
of jumps occurring after finite times, from far away and after long time, so that spatial 
distances and time intervals of normal size appear very small, being x replaced by x/h, t 
replaced by t/i^^/f^ h°'/^) . 
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APPENDIX A: THE CAPUTO TIME-FRACTIONAL DERIVATIVE 

For readers' convenience, here, we present an introduction to the Caputo fractional deriva- 
tive starting from its representation in the Laplace domain and pointing out its difference 
from the standard Riemann-Liouville fractional derivative. In so doing we avoid the sub- 
tleties lying in the inversion of fractional integrals. 

If f{t) is a (sufficiently well-behaved) function with Laplace transform C {f{t);s} = 
f(s) = e f{t) dt , we have 

c[^^f{t);s]=s^f{s)-s^-^f{<d+), 0</3<l, (Al) 



if we define 



We can also write 



dt^-"- ' ' r(l - 13) Jo dT {t - tY 



Eqs. (A. 1-4) can be extended to any non integer /? > 1 , (see e.g. the survey by Gorenflo &l 
Mainardi0) . We refer to the fractional derivative defined by (A. 2) as the Caputo fractional 
derivative, as it was used by Caputo for modelling dissipation effects in linear viscoelasticity 
in the late sixtiesSiUl. 

This definition differs from the usual one named after Riemann and Liouville, given by 
the first term in the R.H.S. of (A. 4), and defined e.g. in the treatise on Fractional Calculus 
by Samko, Kilbas & MarichevL3. 

Gorenflo & MainardS and PodlubnylH have pointed out the usefulness of the Caputo 
fractional derivative in the treatment of differential equations of fractional order for physical 
applications. In fact, in physical problems, the initial conditions are usually expressed in 
terms of a given number of boundary values assumed by the field variable and its derivatives 
of integer order, despite the fact that the governing evolution equation may be a generic 
integro-differential equation and therefore, in particular, a fractional differential equation. 
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APPENDIX B: THE RIESZ SPACE-FRACTIONAL DERIVATIVE 



If f{x) is a (sufficiently well-behaved) function with Fourier transform 

/ + 00 
-oo 

we have 

J" l-fr-f{xy,ti} = -\^rf{K), 0<a<2, {B.l) 



d\x\' 



if we define 



sin^avr/2) ^ /(x + g) - 2/(x) + /(x - Q 
^/(x)=r(l + a) dC. {B.2) 

The fractional derivative defined by (B.2) can be called Riesz fractional derivative^ as it 
is obtained from the inversion of the fractional integral originallv introduced by Marcel 
Riesz, known as the Riesz potentia^^. The representation (B.2)il, is more explicit and 
convenient than others found in the literaturJBil. It is based on a suitable regularization 
of a hyper-singular integral. 

For a = 2, the Riesz derivative reduces to the standard derivative of order 2, as — |kP = 

For a = 1, the Riesz derivative is related to the Hilbert transform, resulting in the formula 



d\x\ vr dx J-oo X — ^ 

We note, by writing — = — (k^)"/^, that the Riesz derivative of order a can be 



interpreted as the opposite of the a/2 power of the (positive definite) operator —D^ — 



namely 

d" ( d 



2 \ "/2 



BA) 

_ rent notation which 1 

into account asymmetries was used by Gorenflo & Mainardi0llll. 



d\x\^ \ dx"^ j 

The notation used above is due to Saichev &; Zaslavsky0. A different notation which takes 
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